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SECTION  I 


INTRODUCTION 


This  report  presents  the  results  of  an  evaluation  of  the  cap¬ 
ability  of  the  short  and  long  period  Norwegian  Seismic  Array  (NORSAR)  to  dis 
criminate  between  Eurasian  earthquakes  and  presumed  explosions.  It  com¬ 
plements  the  analyses  of  the  individual  performance  of  various  discriminants 
presented  in  Special  Report  No.  11,  (Ringdal  and  Whitelaw,  1973)  and  Special 
Report  No.  12,  (Laun,  Shen,  and  Swindell,  1973)  for  the  NORSAR  short  and 
long  period  arrays,  respectively. 

The  NORSAR  array  is  centered  about  109  km  north  of  Oslo, 
Norway,  at  a  latitude  of  60.  8°N  and  a  longitude  of  10.  8°E.  It  consists  of  22 
subarrays,  each  with  6  short  period  and  one  3-component  long  period  sensors 
(Figure  1-1).  The  diameter  of  each  subarray  is  8-10  km;  the  full  array  has 
an  aperture  of  approximately  100  km. 

The  results  presented  in  this  report  are  based  upon  Eurasian 
earthquakes  occurring  between  30  April,  1971  and  31  July,  1972  and  Eurasian 
presumed  explosions  from  1 97  1 ,  1972,  and  1973.  Essentially,  the  data  base 
consists  of  events  that  were  processed  for  both  the  short  period  and  long  per¬ 
iod  evaluations  (Special  Reports  No.  11  and  12,  1973').  However,  a  few  addi¬ 
tional  presumed  explosions  from  1972  and  1973  have  been  analyzed  in  order 
to  increase  the  common  event  population.  The  complete  data  base  is  listed  in 
Tables  1-1  and  1-2,  Total  number  of  events  is  257,  24  of  which  are  presumed 
explosions.  A11  of  the  earthquakes  are  either  shallow  or  of  unknown  depth. 

Figure  1-2  presents  a  breakdown  of  the  processed  events  by  in¬ 
formation  source.  Our  events  have  been  selected  from  seismic  bulletins  pro¬ 
vided  by  four  different  organizations: 


1-1 


o01C  (9) 

A  ° 

O  O  O  , 

O  ' 

A  O 


02 C  (10) 
o 


o  o 

°  03C  (11) 
A  o 


o  o 

o  L  14C  (22) 


o  04 C  (12) 

o  o 


o  O  13C  (21) 
D  * 


01B  <2)0A  Q 


o  o02B  (3) 


o  o  07 B  (8) 

^  o 


O  o  A 


°  O  o 

A  05 C  (13) 


o  D  03B  (4) 


O  A 


o  o  q01A(1) 


o  o  Q12C  (20) 


o  o 


O  o  06B  (7) 
A  ° 


°  06C  (14) 

o 


o 

°  o  ° 

a  lie  (19) 

O 


o  04B  (5) 


o  o 


O  05 B  (6) 


lgC  (18) 


o  o 


o 

0  O  A 

q  A  08C  (16)  °  o  07 C  (15) 


20  km 


0°9C  (17) 


o  Short  Period  Sensor 


a  Short  Period  and  3 -Component  Long  Period  Sensor 


FIGURE  1-1 


NORSAR  SHORT  AND  LONG  PERIOD  ARRAY 


{ 

1 


TABLE  1-1 
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4  7.  °N 

16. 2F 

1  1 

4.0 

p 

TA0/005/12N 

01/05/7? 

1  2.0?. 54 

3  7.  »N 

73.19 

NPR 

4.5 

s 

KOM/005/14N 

01/05/72 

1  4 . 2  f  .  4  P 

56. 6*' 

169. 4F 

MPC 

4.0 

s 

NS  Ml 

KA  v/005  /  16*' 

>1/05/7? 

16.09.50 

5  7 . 7  M 

1 60. 5 F 

9  5 

3.9 

S 

TA  I/006/06N 

01/06/7? 

04.73.34 

2  3.3* 

123.49 

NPR 

4.7 

P 

M  9'N 

SWR/007/20M 

01/07/7? 

? 0.  ?  7. 3? 

44.  IN 

45.19 

NUP 

4.? 

s 

Nl 

KPM/ 00° /03N 

01/09/7? 

03.23.06 

54.4* 

164.49 

NPR 

3.6 

s 

NS  *’l 

K AM/009/14N 

01/09/7? 

14.00.59 

55.7|S 

163.6F 

NPR 

4.7 

s 

KUR/009/14N 

01/09/7? 

I4.47. 46 

45.1* 

14P.4F 

NnK 

7.3 

4 

*S  Ml 

KCM/011/0PN 

01/11/7? 

0P.'4. 

54.7* 

169.29 

7  R 

3  .9 

S 

Nl 

C.P  F/0 12/ 1 3N 

01/12/7? 

17.P1.2C 

7  5.0* 

7  3.5F 

NPR 

4.9 

P 

M  p  NJ 

KAM/012/2OM 

01/12/7? 

20.20.15 

55.  6* 

1 63 . 9F 

NPR 

6.« 

c 

SIH/013/17N 

01/n/7? 

1  7 . 2  4  „  0  7 

61. 9N 

147. IF 

MPO 

5.3 

S 

SI P/014/01N 

01/14/72 

03.2C. i 0 

67.  ^N 

171 .5T 

NPR 

3.9 

D 

VPN 

I R A/014/??M 

01  /  ]  4 / 7  2 

??. 1C. 04 

72.  a*! 

4  6.99 

Njr.  u 

6.1 

D 

KIIP/01  5/OON 

01/16/7? 

00. 5P. 37 

4  9.6  *' 

155. 09 

NpR 

7.9 

s 

NS  Ml. 
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FVFNT 

1  M  f  1  N 

s  4"  if  r  5 

OES  IGNAT  I  ON 

0  A  T  F 

T  I  Mr 

1  A  T 

1  ON 

HFPTH 

yc 

t>l  T  \ 

r .  'MM TMT 

IP  A/0IP/21N 

01/1F/7? 

21.12.C2 

3  7 . 4  N 

4  0 . 7  5 

\r'P 

4.4! 

P 

M  fl\: 

nnn/o?o/o?N 

01/20/7? 

0?.  15.07 

36.  sV 

>7.  IF 

NH7 

4.  o 

s 

KHP/02P/01N 

01/22/72 

01.41.24 

C0.0M 

15?. 04 

N <r  o 

4.? 

p 

1  n\j 

TIIR/022/17M 

01/22/7? 

17.17. 

0  7.  f-l 

?C.  QP 

1  ? 

4  .  4 

n 

KAM/025/10N 

01/25/7? 

1 0.02. 40 

5  ?  .  5  * 

1  50.  Of 

N?P 

4  .* 

D 

M| 

GRF/026/1 ?N 

01/25/7? 

1  2 . 64 • ?° 

?4.  5N 

?  5 .  5  F 

NPP 

4.0 

9 

KAM/027/20N 

01/27/72 

?0.?7. ?P 

55.  7 N 

16?. 35 

40 

J  .5 

5 

VI  | 

FCS/028/O4N 

01/?°/ 7? 

C  4  .  ?  ?  .  ?  o 

2  7.5  V 

126.65 

V  n  K 

4.4 

F 

K IP/028/20N 

01/28/72 

.  0 . 2  c .  1  ° 

4?. ON 

7  P  •  0  5 

NOP 

4.4 

P 

M  RM 

FRS/028/21N 

91/25/72 

2 1 . qc.or 

4  5.0V 

1  ?6.  Q5 

vr)7 

4.0 

V 

V'l 

KUP/028/23N 

01 /?«/?? 

2’. *2. 51 

4  ° .  ?  N 

1  57.  K 

mop 

3.° 

6 

IP  A/029/0°N 

01/29/7? 

C°. tC. 5  P 

2  7.  ON 

62.05 

N‘|'0 

3.0 

n 

KAV/032/10N 

02/01 / 72 

10.15.C° 

5  5.  H  V 

1 5?. PF 

NPc 

4.1 

6 

ML 

KUR/033/09N 

02/02/7? 

0°. q P. 51 

4  6.  3  N 

1 4  6  o  4  F 

N”F 

?  .6 

6 

V'l 

KIIR/0’3/  17N 

9  2/02/7.? 

1  7 . 5  5  .  3  r 

5  0.  7  N 

1 6,  o .  i  r 

NOP 

?.6 

F 

ML 

YIJN/034/07N 

■12/CV77 

0?. ??. 40 

2  'AN 

1C?. 45 

vtrr 

4.5 

P 

RAI/O?.  5/03N 

02/04/7? 

0*. 34. 54 

5  1 . 4  N 

1  1  5  •  0  r 

MOF 

4.2 

C 

IT  A/035/0 9N 

0?/04/7? 

0°. 1 p. V 

4?  AN 

1  ?• .  2  F 

?  3 

4.4 

P 

IRA/041/09N 

02/10/7? 

0c.C4.0O 

29. 6N 

6  P •  OF 

NT  5 

?  .0 

0 

I» A/041/16M 

92/10/7? 

1  5.40.  15 

?  9 . 5  N 

5  0  .  °  5 

tl 

4.1 

p 

S  IVI/042/09N 

o;/n  /  7  ? 

0  5  4  / 

?°.°N 

7  7.45 

*>  "3 

4.4 

p 

T IP/042/1 2M 

02/11/7? 

1  2 • 2 1 • 4° 

2  5.  ON 

“  7.0° 

NF’P 

4.7 

Vi 

V  1 

KAM/0A2/21N 

02/11/7? 

?1  .  ?5  .  1  7 

56.1  N 

167. cr 

44 

4 .6 

P 

KUR/044/05N 

o?/n/7? 

05.24.57 

4  ?.  r>  N 

1 47. OF 

NT-  p 

T  Cl 

• 

6 

Ml 

GPF/044/13N 

93/1 ?/7? 

1’ .07,11 

?7.  1  N 

?  4 . 0r 

7  7 

4.5 

P 

KC1M/G44/22N 

02/1  V72 

??.  V  .54 

55.  ?N 

HA-.  6  C 

VJOP 

3.° 

c, 

KIJP/04fc/16N 

02/15/7? 

1  /  .  4  5 . 2  ? 

4  5.0V 

157.05 

NF-i? 

4.1 

5 

GPF/047/OON 

02/15/7? 

00.42.24 

3  6.5  v 

24.  ?  F 

NMF 

4.6 

P 

KUR/049/18N 

0  2  /  1  R  /  7  2 

1°.C?.?4 

4  ?  .  6  N 

147. a- 

36 

4.7 

5 

SIN/051/1  ON 

0?/?n/72 

10.22.45 

7  «  .  5  N 

°C.  55 

1  6 

3.° 

1 

■ICJf.l 

KAP/051/2CN 

02/20/7? 

20.06.11 

50.  PV 

1  **  1 . 5  r 

MfiP 

4  .  1 

! 

KAM/052/22N 

02/21/7? 

2?  .00. rQ 

5  4 .  4  V1 

161.75 

n  n  r 

4 

I 

YU G/ 052/2 3 N 

02/21 /7? 

?  1 . 0  r .  q  c 

4  1 . 0  N 

22.35 

VP  Q 

4.0 

i 

V  1 

M0N/053/01N 

02/22/72 

0  1  .*-?.?  ( 

4°.  ON 

11  5. 

V  "  P 

4.1 

I 

'  l 

HIN/053/0PN 

02/22/7? 

CP. 1 4. ?6 

3  6.6  N 

6  p  .  5  5 

NPF 

4  .0 

1 

VPM 

KIJR/054/03N 

0  7  / ,?  3  /  7  ? 

f  ?  .  4  ?  .  4  1 

4  3.4V 

1  4  F'  .  7  6 

7  0 

4.0 

I 

K  A  M/054 / 1 9N 

9?/  ? / 7? 

1  c  .3  7. 20 

5  5  .  ()  N 

1 63. nr 

Nrr 

3.7 

I 

MC 

V'l 

KUP/055/10N 

02/24/7? 

10.1 °. 77 

4  0  .  R  N 

155. 7 c 

M  i  >  u> 

5.4 

1 

K UR/055/ 1 PN 

0?/?5/7? 

15.17. 34 

4  C.'ll 

15  -'.04 

NPF' 

3.5 

1 

"F 

V  1 

KUP/056/19N 

02/?t'/7? 

1  4.  2° 

46.0V 

14  7.04 

N5P 

^  c 

• 

T 

r  i 

WPS/056/22N 

0?/?q/72 

2  2  .  ?  4 .  4  0 

50.  ON 

?".  Or 

fqr 

0  *7 

•  ■ 

1 

r  l 

KUR/056/22N 

92/2q/7^ 

22.43.07 

4  6.7V 

156 . 0r 

N45 

4  .41 

1 

KUR/057/05N 

02/26/7? 

CK.r-F.  2? 

46.  p  N 

15?. 65 

ivnc 

4.4 

I 

KAM/057/09M 

02/25/7? 

0C.0* . ’? 

55.0V 

1 1”  2.  OF 

V  HP 

?.? 

1 

*•  5 

f  1 

ERS/057/  15N 

02/25/7? 

15.06.4? 

5  9 . 3  V 

l?  5.7T 

vpr 

7  .« 

T 

•  P\j 

Ml 

1-5 
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FVENT 

OR  IGIN 

sou 

POE 

DES  IGNAT  ION 

DATE 

TIME 

LAT 

LON 

OF  PTH 

MR  RL 

TN  COMMENT 

YUN/057/18N 

02/26/72 

18.56.13 

27. IN 

100. OF 

NDR 

4.7  I 

L01/058/08N 

02/27/72 

08.4?. 57 

88. ON 

74.  OW 

NOR 

3 .3  I 

NS 

Ml 

LOM/05  8/ 1CN 

02/27/72 

10.03.03 

87. ON 

53.51* 

NOR 

4.9  I 

•  L 

L0M/058/11N 

02/27/72 

11.03.19 

90.  ON 

95.01* 

NOP 

3.5  I 

M| 

LCM/058/17N 

02/27/72 

17.50.25 

86. 2N 

77. 2W 

NOR 

4.4  I 

BAI/058/22N 

02/27/72 

22.15.03 

55. ON 

93. 2F 

NOR 

4.5  I 

KUR/059/01N 

02/28/72 

01.04.22 

46.  ON 

148. OE 

NOR 

4.2  I 

NS 

NL 

PAK/059/05N 

02/28/72 

C8.1P.56 

36.  7N 

71. 4E 

NOR 

4.2  I 

MRN 

KAM/059/ I IN 

02/28/72 

11.35.31 

56.  ON 

163. OF 

NOR 

4.1  T 

NS 

ML 

AFG/059/ 18N 

02/28/72 

18.12.35 

36. ON 

68. 7F 

NOR 

4.4  I 

MRM 

KAM/059/ 20N 

02/28/72 

20.04.00 

56. IN 

164. 2E 

NOR 

3.6  I 

NS 

NL 

I RQ/060/08N 

02/29/72 

08.02.51 

32. 8N 

46. 6E 

NOR 

4.0  I 

MRN 

Mi  L 

IR  A/062/1  AN 

03/02/72 

14.1C.13 

31. 6N 

42. IE 

NOP 

4.0  I 

MRN 

NL 

ALM/062/19N 

03/02/72 

19.57.42 

43. ON 

76. OF 

NOR 

3.5  I 

MRN 

NL 

KAM/063/00N 

03/03/72 

00.39.23 

53. ON 

159. ?E 

NOP 

4.1  I 

NL 

K0M/063/0PN 

03/03/72 

08.13.55 

55.  8N 

163. 9E 

NOP 

4.1  I 

NL 

YUG/063/2IN 

03/03/72 

21.26.51 

44.  7N 

1  8.4E 

3? 

4.4  I 

ML 

KUR/063/23N 

03/03/72 

23.10.41 

50. 2N 

155. 7E 

NOR 

4.5  I 

SIN/ 064/04N 

03/04/72 

04.00.09 

40. 2N 

79. OE 

NOR 

4.5  I 

KAM/066/0AN 

03/06/72 

06.05.08 

53. 5N 

160. OE 

NO  « 

3.9  I 

NS 

NL 

KUR/066/09N 

03/06/72 

09.55.09 

45. ON 

150. OF 

NOR 

3.7  I 

NL 

0KH/066/I9N 

03/06/72 

19.13.25 

56. ON 

140. OE 

NOP 

4.2  I 

ML 

CHI/066/23N 

03/06/72 

23.17.53 

40. ON 

103. OF 

NOR 

4.5  I 

NS 

L. 

NL 

YUG/067/05N 

03/07/72 

05.21.21 

43. ON 

21.  OF 

NOP 

2.7  I 

MRN 

NL 

OKH/068/02N 

03/08/72 

02.38.11 

51 ,2N 

151. 9F 

NOR 

4.2  I 

NL 

IRA/068/2IN 

03/08/72 

2i.45.il 

27. 6N 

56.  7E 

45 

4.6  I 

PUL/068/22N 

03/08/7? 

22.04.0 2 

40.  RN 

22. 8E 

NOB 

3.5  I 

MRN 

KUR/070/06N 

03/10/7? 

06.50.18 

45.  IN 

149. 5E 

NOR 

3.7  I 

ARC/ 071/06N 

03/11/72 

06.47.07 

81. ON 

157. OE 

NOR 

4.3  1 

NS 

NL 

KAS/07 1/ 13N 

03/11/72 

13.31.39 

35. ON 

76. OE 

NOR 

4.1  I 

MRN 

Nl 

KUR/073/02N 

03/13/72 

02.11.05 

49.  ON 

158. OF 

NOR 

3.8  I 

NL 

AFG/07  3/05N 

03/13/7? 

05.45.13 

37. ON 

70. OF 

NOP 

4.0  I 

MRN 

NL 

TIB/073/18N 

03/13/72 

1°. 27.07 

34. ON 

83. OE 

NOR 

4.1  I 

MRN 

NL 

TIB/075/06N 

03/15/72 

06.00. 33 

30. 4N 

84. 5E 

NOP 

5 .3  I 

MRN 

KUR/077/07N 

03/17/72 

07.48.02 

49. ON 

156. 2E 

NOR 

5.2  I 

TAC/077/09N 

03/17/72 

08.17.  11 

40. IN 

69.  7E 

26 

5.2  T 

IRA/077/I7N 

03/17/72 

17.1!. 28 

28. ON 

84.  OE 

NOR 

3.9  ! 

MRN 

NL 

KAS/077/23N 

0V17/72 

23.33. 37 

32. ON 

75. OE 

NOR 

3.5  1 

MBN 

NL 

KAZ/078/07N 

03/18/7? 

07. 11. 55 

47. ON 

81. OE 

NOR 

3.6  I 

MRN 

K AM/078 / 1 3N 

03/18/7? 

13.52.14 

57. ON 

l 61 . OF 

NOP 

3.6  T 

NL 

KAM/078/18N 

03/18/72 

18.?c.37 

50. 6N 

156. 7E 

Nnp 

4.7  I 

OKH/078/I9N 

03/18/7? 

19.17.25 

54. ON 

150. OE 

NOP 

3.7  I 

NL 

CAU/079/03N 

03/19/7? 

03.34.31 

42. 7N 

38. IF 

NOR 

3 .9  I 

KUR/080/14N 

03/20/7? 

14.08.12 

47. ON 

154. OE 

NOP 

4.0  I 

MBN 

NL 

SIN/080/2IN 

03/2C/7? 

21.47.55 

40. ON 

80. OE 

NOR 

3.4  I 

MRN 

nl 

*  a- 


EVFNT 

OEblGNAT ION 

CKZ/153/01N 

MON/  1  *53/  1  IN 

S IN/154/05N 

SIN/154/06N 

T  S I / 1 54/ 1 6N 

CHI/ 154/ 20N 

RYU/155/02N 

IRA/I55/08N 

IRA/156/03N 

KAM/I56/07N 

K  A  M/  I  57/04N 

GR  E/ 157/1  ON 

IRA/ 157/1  IN 

P AK/ 157/ 1 IN 

NFJ/157/19N 

KIJR/15B/06N 

TAI/160/10N 

TUR/160/12N 

RUR/160/16N 

TIA/160/23N 

CRE/161/07N 

P AK / 162/ 1  IN 

KUI/163/23N 

KII2/163/23N 

KUR/  164/00N 

I P 1/ 1 66/ 13N 

IRA/ 165/OON 

KAM/165/04N 

CAS/166/00N 

IRA/166/04N 

IR2/166/12N 

GRE/167/00N 

KAM/168/09N 

KAM/ 168/22N 

IP A/16R/23N 

AUS/169/09N 

KUR/ 169/1 9N 

TIR/170/04N 

KUR/i 70/09N 

KIJ  2/  171/1  8N 

KUR/171/22N 

ERS/172/09N 

KAS/172/15N 

TIJR/173/05N 

KAM/ 17  3/ ION 
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DATF 

T  IMF 

LAT 

1  ok 

Pc  nTH 

NR 

Ml  TN 

rOMNFK  T 

06/01/72 

01 .??. 26 

52.  ON 

70. OF 

KOR 

7  .6 

S 

NS 

NL 

06/01/7? 

11.??.  15 

44.  OK 

1 03. or 

Nnp 

3.7 

N 

06/02/72 

C5.ll . 13 

43. ON 

Pl.Of 

NOR 

3.5 

N 

Ml. 

06/C2/7? 

C  6  .  ?  C  • 49 

42. OK 

0  l  .  OF 

KH« 

7.0 

K 

06  / 0 2 / 7 2 

16.49.2? 

36. OK 

92.  or 

NOR 

7.7 

N 

06/C2/72 

?0.3?. 55 

2  8.4N 

c5.of 

NOP 

4.3 

o 

06/03/72 

02.16.51 

23.  SN 

125. 5F 

NOP 

5.0 

P 

06/03/72 

08.21.30 

25.  OK 

5  3 . 0  p 

NOR 

4.? 

K 

06/04/72 

07.37. 45 

3C.OK 

*4.  OF 

NOR 

4.2 

N 

06/04/72 

07.62.38 

53. ON 

158.  or 

NOP 

4.0 

S 

06/05/72 

04.12.54 

5  6.  2  N 

168. if 

NT  R 

4.3 

P 

06/05/7? 

10.44. 5C 

3  7.  7N 

?  1 . 4F 

60 

4.7 

P 

06/05/72 

11.17.57 

34.  ON 

46.  OF 

NOR 

3.9 

K 

06/05/72 

11 . 52.  c3 

2  5.  8N 

70.  3F 

27 

4  .8 

P 

06/05/7? 

19.CC. 1? 

86.  5 K 

"’P.OF 

NOP 

4.k' 

5 

06/06/72 

06.32.  10 

40.  ON 

155. OF 

NOP 

7  .6 

s 

NS 

Nl 

06/08/72 

10. 17.44 

21. OK 

120.29 

NOR 

4.9 

p 

06/C6/7? 

12.46.15 

41. ON 

44. OF 

NOR 

4.1 

s 

06/08/72 

16. CP. C6 

19. OK 

r4.  OF 

NOP 

4.3 

M 

06/08/7? 

23. 1C. 1? 

29. 5n 

o?.3F 

64 

4.7 

p 

06/09/7? 

07.42.20 

34. 8K 

26. 5F 

KOP 

4.0 

p 

06/1C/7? 

11.29.11 

28. 2K 

66. 5F 

KPR 

4.5 

p 

06/11/7? 

?”1.2'J.C4 

48.  ON 

15?. OF 

NOR 

4.0 

s 

Ml. 

06/11/72 

23.33.44 

47. ON 

1 c  ? • OF 

NOP 

4.3 

5 

Nl 

06/12/72 

00.15.16 

44. ON 

14P.0P 

NOR 

3.7 

5 

Nt 

06/12/72 

13.34.01 

3  3.1K 

46. 3F 

NOR 

5.4 

P 

06/13/72 

00. 5C. 37 

33.  IN 

46. 3F 

27 

5.1 

P 

06/17/7? 

04.53. 30 

55.  OK 

162. OF 

NOR 

7.8 

S 

Ml 

06/14/7? 

00.45. 54 

40.  IN 

51. OF 

4  7 

4.7 

P 

06/14/7? 

04.34.28 

3  3.  ON 

46. IF 

KOP 

5.7 

p 

06/14/72 

12.35.C5 

2  7 . 0  K 

56. 0F 

NOR 

7.6 

N 

NL 

06/15/7? 

00.33.24 

38. 3N 

?  2 . 2  r 

7  6 

4.0 

P 

06/16/72 

09.54. 41 
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ABBREVIATIONS 

P  :  Reported  by  PDE 

I  :  Reported  by  ISM 

S  :  Reported  by  SDAC/  LASA 

N  :  Reported  by  NTNF/NORSAR 

NS  :  No  short  period  signal  detected  on  NORSAR  data 
(TI  analyst  decision) 

NL  :  No  long  period  signal  detected  on  NORSAR  data 
(TI  analyst  decision) 

MBN  :  Bodywave  magnitude  is  based  on  NORSAR  data 

NOR:  Normal  depth  -  means  that  focal  depth  has  been  restricted 

at  33  km;  i.  e.  there  is  not  sufficient  evidence  (such  as  depth 
phases)  to  provide  a  specific  depth  estimate. 


•  The  PDE  listings  (Preliminary  Determination  of  Epicenters) 
issued  by  the  National  Oceanic  and  Atmospheric  Administra¬ 
tion  (NOAA) 

•  The  LAS  A  seismic  bulletin  issued  by  the  Seismic  Data  Analysis 
Center  (SDAC) 

•  The  NORSAR  seismic  bulletin  compiled  at  Kjeller,  Norway 

s  The  bulletin  from  the  International  Seismic  Month  (ISM),  which 

covers  February  2  0  to  March  19,  1972,  and  has  been  compiled 
at  Lincoln  Laboratories. 

In  all  cases  where  an  event  was  repcrted  by  more  than  organi¬ 
zation,  we  selected  our  source  information  according  to  the  priority  list:  ISM, 
PDE,  LASA,  NORSAR.  As  can  be  seen  from  Figure  1-2,  most  events  of  m 

b 

above  4.5  have  been  reported  by  PDE,  while  LASA  and  NORSAR  "upply  the 
detection  information  for  most  low  magnitude  events. 

Figure  1-3  is  a  histogram  showing  the  distribution  of  processed 
events  as  a  function  of  epicentral  distance  from  NORSAR,  and  indicating  which 
events  are  presumed  explosions.  It  is  seen  that  almost  all  our  presumed  ex¬ 
plosions  occur  at  a  distance  of  about  38  degrees,  this  corresponds  to  the 
Eastern  Kazakh  region  in  the  Soviet  Union.  With  respect  to  close-in  pre¬ 
sumed  explosions  from  Western  Russia,  the  number  of  events  included  in  this 
report  is  much  lower  than  the  corresponding  number  for  our  SP  evaluation. 

The  reason  for  not  processing  more  of  these  events  is  mainly  the  occurrence 
of  amplitude  clipping  for  these  large  signals,  which  has  made  any  measure¬ 
ment  of  surface -wave  magnitudes  impossible  in  many  cases. 

The  earthquake  population  in  our  data  base  may  conveniently 
be  divided  into  three  regions  where  the  majority  of  events  are  located: 


•  The  Mediterranean  area  (Distance  15-30  degrees).  Approxi¬ 
mately  30  events  from  this  region  have  been  processed. 

•  Iran  and  Central  Asia  (Distance  35-50  degrees).  The  number 
of  events  is  approximately  70. 

•  Kamchatka  and  the  Kuriles  Islands  (Distance  60-70  degrees). 
Approximately  100  events  from  this  reg'on  are  included  in  our 
data  base. 

A  total  of  nine  discriminants  have  been  evaluated  in  this  report. 
Five  of  these  are  based  on  short  period  data  only;  the  remaining  four  combine 
short  and  long  period  information.  The  discriminants  are  briefly  described 
as  follows: 

1.  P30  Mean  Square 

This  discriminant,  which  is  a  measure  of  event  complexity,  is 
computed  by  crosscorrelating  4  seconds  of  the  waveform  (beginning  a  few  points 
before  P-wave  onset)  with  the  next  30  seconds  of  the  waveform  and  with  the 
noise  preceding  the  signal.  A  mean  square,  weighted  by  the  lag,  is  then  com¬ 
puted  from  the  correlations  over  both  30  seconds  of  the  noise  and  30  seconds 
of  the  signal.  The  noise  mean  square  is  subtracted  from  the  signal  mean  square 
to  obtain  the  discriminant  used  (Texas  Instruments  Incorporated,  1971). 

2.  Autocorrelation  Mean  Square 

This  discriminant  is  also  a  measure  of  complexity.  The  auto¬ 
correlations  of  a  30-second  noise  gate  and  of  a  30-second  signal  gate  are  com¬ 
puted  and  a  weighted  mean  square  then  derived  from  these  correlations  for  the 
noise  and  signal.  The  discriminant  is  derived  from  the  signal  mean  square 
minus  the  noise  mean. 


Envelope  Difference 
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3. 

This  discriminant  is  also  derived  from  tl  r  P30  correlation  by 
computing  the  mean- square  difference  between  the  envelope  correlation  and  a 
fixed  decaying  exponential,  the  decay  rate  of  which  is  the  average  rate  for  an 
ensemble  of  16  explosions  recorded  at  LASA.  As  with  the  first  two  statistics, 
envelope  difference  is  a  measure  of  complexity. 

4.  Dominant  Period 

This  discriminant  is  computed  by  finding  the  cycle  in  the  wave¬ 
form  with  the  maximum  absolute  amplitude;  the  dominant  period  is  the  duration 
of  this  cycle  in  seconds.  This  parameter  can  be  estimated  with  some  confidence 
even  for  events  with  a  relatively  low  signal-to-noise  ratio.  The  dominant  period 
discriminant  is  a  rough  measure  of  spectral  energy  distribution. 

5.  Spectral  Ratio 

This  discriminant  is  derived  from  the  signal  power  spectrum 
over  a  gate  beginning  just  before  the  signal  arrival.  The  power  spectrum  is 
smoothed  over  three  frequency  points,  and  the  power  in  three  bands  is  com¬ 
puted;  Band  1:0-0.  55  Hz;  Band  2:  0.  55  -  1.  5  Hz;  Band  3:  1.  5  -  5.  0  Hz. 

These  bands  have  been  selected  based  on  NORSAR  data.  Spectral  ratios  com¬ 
puted  were  Band  3  to  Band  2  and  Band  3  to  Band  1,  respectively. 

6.  (Rayleigh)  -  Discriminant 

This  discriminant  compares  the  magnitude  of  the  Rayleigh 

wave  signal  measured  at  NORSAR  to  the  reported  bodywave  magnitude  of  each 

event.  The  detailed  procedure  of  measuring  M  is  described  in  Special  Re- 

s 

port  No.  12  (1973). 

7.  M  (Love)  -  m,  Discriminant 

s  b 

This  discriminant  compares  the  magnitude  of  the  Love  wave 
signal  measured  at  NORSAR  to  the  reported  bodywave  magnitude  of  each  event. 


8. 


AR/m  Discriminant 
b 


The  AR  value  is  a  measure  of  the  Rayleigh  wave  energy  of  an 
event  (Brune,  Espinosa  and  Oliver  (1963);  Evernden  (1969)).  Our  AR  values 
have  been  normalized  to  a  bodywave  magnitude  of  5.  0  as  described  in  Special 
Report  No.  12  (1973). 


AL/ rn  Discriminant 
b 


The  AL  values  are  computed  by  the  same  method  as  described 
for  AR,  except  that  in  this  case,  measurements  are  taken  of  the  Love  wave 


energy. 


All  the  nine  discriminants  defined  above  are  two-dimensional. 


Ifi  the  sense  that  discriminant  values  are  plotted  against  bodywave  magnitudes 


(mb)  in  order  to  obtain  a  good  separation  between  earthquakes  and  presumed 


explosions.  It  is  thus  important  to  obtain  an  m,  estimate  which  is  as  reliable 

b 


as  possible  for  each  event.  We  have  chosen  to  use  the  value  quoted  in  the 


source  bulletin  (PDE,  ISM,  LASA,  NORSAR)  in  almost  ail  _ases;  the  only  ex¬ 


ception  is  events  for  which  the  PDE  mfa  has  been  based  upon  near -field  stations 


only.  In  those  cases,  we  have  used  the  NORSAR  m,  values. 

b 


The  first  five  discriminants  defined  above  were  applied  to  all 
events  with  SP  detections  at  NORSAR;  similarly,  the  last  four  were  evaluated 
for  all  LP  cetections.  Table  1-3  summarizes  the  number  of  detected  events; 


it  is  seen  that  about  90  percent  of  all  reference  events  were  detected  on  SP 
data,  while  about  70  percent  had  LP  detections  (Rayleigh  waves,  Love  waves 
or  both).  The  NORSAR  LP  detection  capability  is  actually  somewhat  lower 
tnan  these  numbers  indicate,  since  "mixed"  events  (i.  e.  ,  events  for  which  the 
LP  signals  have  been  buried  in  signals  from  a  stronger  event)  were  not  includ¬ 
ed  in  our  data  base.  It  was  found  in  Special  Report  No.  12  (1973)  that  about 
20  percent  of  the  events  originally  selected  for  that  report  were  of  this  type. 
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DETECTION  PERFORMANCE  OF  THE  NORSAR  SP  AND  LP 
ARRAYS  FOR  THE  REFERENCE  EVENT  SET 


SECTION  II 

NORSAR  DISCRIMINATION  RESULTS 


A.  INTRODUCTION 

This  section  presents  the  results  achieved  by  applying  the  nine 
discriminants  defined  in  Section  I  to  NORSAR  data.  The  performance  of  each 
individual  criterion  is  described  in  Subsection  II-B,  while  combined  criteria 
are  evaluated  in  Subsection  II-C. 

\  method  for  measuring  the  separation  achieved  by  single  or 
multiple  discriminants  '/as  introduced  in  Special  Report  No.  11,  1973,  and 
will  also  be  applied  here.  Basically,  the  method  consists  of  finding  the  best 
linear  separation  between  two  point  sets  in  a  plane,  in  the  following  way 
(Figure  II- 1): 

1.  For  any  g:ven  straight  line,  the  distances  from  the  line  are 
computed  for  all  points  corresponding  to  events  in  the  earthquake  and  pre¬ 
sumed  explosion  populations. 

2.  The  two  sets  of  real  numbers  thus  obtained  are  considered  as 
sampled  from  two  Gaussian  populations;  ore  N(fij,  (7^)  (earthquakes)  and  one 
N(J1_,  (T  )  (presumed  explosions). 

£t  L 

3.  The  straight  line  that  gives  two  Gaussian  populations  with  the 
best  separation  is  chosen  and  the  measure  of  separation  is  defined  as  the 
corresponding  probability  of  correct  classification. 

There  is  no  unique  method  of  computing  the  separation  of  two 
Gaussian  distributions  with  unequal  variances.  We  chose  under  3.  above  to 
measure  the  separation  by  computing  the  worst  case  probability  of 


Discriminant  Value 


FIGURE  II- 1 

ILLUSTRATION  OF  A  METHOD  TO  MEASURE  THE  SEPARATION 
BETWEEN  TWO  POINT  SETS  IN  A  PLANE 


misclassific ation,  applying  the  well-known  minimax  principle  (van  Trees, 
1968).  The  minimax  separation  threshold  between  the  two  populations  is 
given  as 


T  = 


M2ai  +M1CT2 
°l+<,2 


and  the  associated  probability  of  error  is 

P(miss)  =  P(false  alarm)  = 


( im  \ 

V  °1  +  °2  I 


where  $  is  the  standard  cumulative  Gaussian  distribution  function. 

Thus  it  is  seen  that  for  the  above  threshold  T  ,  the  condition¬ 
al  probabilities  of  correctly  identifying  an  earthquake  and  correctly  identify¬ 
ing  a  presumed  explosion  are  identical,  and  therefore  provide  a  well-defined 
measure  of  the  separation  of  the  two  populations. 

The  method  of  measuring  separation  described  above  may  be 
easily  adapted  to  multivariate  discriminants.  In  fact,  since  the  method  re¬ 
duces  a  two-dimensional  discriminant  to  two  sets  of  real  numbers,  it  can  be 
applied  successively  to  reduce  any  given  array  of  discriminant  values  to  the 
one -dimensional  case. 


B.  PERFORMANCE  OF  INDIVIDUAL  DISCRIMINANTS 

Individual  discriminant  values  for  each  of  the  discriminants 
defined  in  Section  I  are  plotted  as  a  function  of  bodywave  magnitude  in  Figures 
II-2  through  II- 11.  Earthquakes  are  represented  by  a  cross,  while  presumed 
explosions  are  indicated  by  an  asterisk.  All  short  period  discriminants  have 
been  evaluated  on  the  adjusted-delay  array  beam.  The  complete  set  of  dis¬ 
criminant  values  corresponding  to  these  plots  are  listed  in  Appendix  A. 

The  "best  separation  line"  as  defined  In  Subsection  II- A  be¬ 
tween  the  earthquake  and  presumed  explosion  populations  has  been  drawn  for 
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each  discriminant  plot.  Note  that  this  line  has  been  found  on  the  basis  of  all 

events  with  m  >  4.  4  ;  this  is  because  all  our  presumed  explosions  are  of  at 
b 

least  this  magnitude.  For  the  AR  and  AL,  discriminants  (Figures  II- 10  and 
11-11),  the  separation  line  has  been  restricted  to  be  parallell  with  the  m^ 
axis,  since  these  two  discriminants  are  already  normalized  with  respect  to 
bodywave  magnitude. 

For  each  discriminant,  the  distance  from  the  best  separation 
line  is  a  measure  of  how  well  this  discriminant  has  been  able  to  classify  a 
given  event.  These  distances  are  listed  in  Tables  II- 1  and  II- 2  for  presumed 
explosions  and  earthquakes,  respectively.  In  order  to  make  the  numbers  com¬ 
patible  between  discriminants,  all  distances  were  normalized  by  dividing  by 
the  standard  deviation  in  the  respective  earthquake  population.  A  positive 
distance  means  that  the  event  lies  on  the  "correct"  side  of  the  separation  line, 
while  a  negative  distance  indicates  misclas sification.  Figure  11-12  illustrates 
how  this  distance  relates  to  the  probability  distribution  of  the  earthquake  popu¬ 
lation,  which  is  assumed  to  be  Gaussian  In  fact,  the  probability  of 

an  earthquake  exceeding  the  distance  D  from  the  separation  threshold  is 

P(exceed  D)  =  1  - 

where  <t>  is  the  standard  cumulative  Gaussian  distribution  function. 

It  is  seen  from  the  discriminant  plots  and  confirmed  by  Tables 
II- 1  and  II- 2  that  the  four  combined  SP-LP  discriminants  achieve  essentially 
complete  reparation  between  the  earthquake  and  presumed  explosion  popula¬ 
tions.  In  contrast,  the  short  period  discriminants  do  not  perform  nearly  as 
well,  and  all  of  them  misclassify  several  earthquakes  and  several  presumed 
explosions. 
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Distribution  of  distances  from  separation 
line;  earthquake  population 


I 


FIGURE  11-12 

ILLUSTRATION  OF  THE  RELATIONSHIP  BETWEEN 
THE  EARTHQUAKE  PROBABILITY  DISTRIBUTION 
AND  THE  NORMALIZED  DISTANCE  FROM  THE 
SEPARATION  LINE 


Among  the  presumed  explosions,  the  following  events  were 
found  to  be  difficult  to  classify  properly: 

•  KAZ / 356/ 06N;  m  =  6.  0,  M  (Rayleigh)  =4.5,  M  (Love)  =  3.  3 

D  S  S 

Distance  =  25  degrees 

This  presumed  explosion  from  Western  Kazakh  had  an 
unusually  strong  Rayleigh  wave  signal,  and  was  a  borderline 
decision  for  the  two  discriminants  based  on  Rayleigh  waves. 
Long  and  short  period  signal  traces  for  this  event  are  shown 
in  Figures  11-13  and  11-14,  respectively.  The  Love  wave  dis¬ 
criminants  and  the  SP  criteria  worked  well  for  this  event. 

•  WKZ/233/02N;  m^=  5.  7,  M^(Rayleigh)  =  3.  1,  Ms(Love)=2.8 
Distance  =  24  degrees 

This  presumed  explosion  had  an  epicenter  very  close  to 
that  of  the  event  discussed  above,  KAZ/ 356/ 06N;  the  location 
difference  being  only  about  two  degrees.  Long  and  short  per¬ 
iod  signal  traces  for  WKZ/233/02N  are  displayed  in  Figures 
11-15  and  11-16.  There  are  marked  differences  between  this 
event  and  KAZ/356/06N,  in  particular  with  respect  to  the  rela¬ 
tive  amplitudes  of  the  Love  and  Rayleigh  waves.  Also,  when 
comparing  the  short  period  signal  traces  for  the  two  events,  it 
appears  that  WKZ/233/ 02N  has  less  high  frequency  energy  and 
also  has  a  less  impulsive  signal  onset  than  KAZ/356/06N.  In 
fact,  all  the  short  period  discriminants  fail  to  classify 
WKZ/233/ 02N  properly,  while  all  four  SP-LP  criteria  worked 
well  for  this  event. 

•  EK1/ 345/ 04N;  m  =  5.  7,  M  (Rayleigh)  =4.0,  M  (Love)  =  4.  2 
Distance  =  38  degrees 

This  presumed  explosion  was  misclassified  by  four  dis¬ 
criminants.  The  obvious  reason  for  these  problems  is  a  second 


FIGURE  11-14 

SHORT  PERIOD  TIME  DOMAIN  PLOT  FOR  EVENT 
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FIGURE  11-15 

LONG  PERIOD  TIME  DOMAIN  PLOT  FOR  EVENT 

WKZ/233/02N 
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FIGURE  11-16 

SHORT  PERIOD  TIME  DOMAIN  PLOT  FOR  EVENT 

WKZ/233/02N 
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presumed  explosion  from  approximately  the  same  location 
occurring  just  a  few  seconds  later,  with  an  amplitude  about 
twice  as  large  (Figure  11-17).  This  of  course  distorts  the 
complexity  criteria  (except  for  the  autocorrelation)  and  also 
affects  the  M^-m^  relationships.  However,  the  SP  discrimin¬ 
ants  based  on  spectral  content  were  able  to  classify  this  event 
properly. 

•  EKZ/277/01N;  m^  5.3,  M^(Rayleigh)  =  3.5,  Distance  =  38 
degrees 

This  presumed  explosion  was  misclassified  by  three 
short  period  discriminants.  Figure  11-18  is  a  plot  of  the  SP 
signal  traces  for  this  event.  There  is  noticeably  less  high 
frequency  energy  in  this  signal  than  has  been  observed  for 
other  Eastern  Kazakh  events,  e.  g.  ,  EK1/345/04N,  Figure  11-17. 

•  URA/191/16N,  m  =  5.  3;  EKZ/246/08N,  m  =  5.  1;  EKZ/349/07N 

d  b 

"V  4- 9 

No  surface  wave  was  detected  at  NORSAR  for  any  of  these 
three  presumed  explosions.  However,  the  SP  discriminants 
were  able  to  place  all  three  events  in  the  proper  category. 

•  EKZ/188/01N;  m  =  4.4 

b 

This  presumed  explosion  was  misclassified  by  all  SP  dis¬ 
criminants.  Surface  waves  were  not  detected  for  this  event.  It 
is  thought  that  our  SP  criteria  are  not  very  effective  for  tele- 
seismic  events  of  a  magnitude  this  low. 

With  respect  to  earthquakes,  we  proceed  to  briefly  discuss 
those  above  mb  =  4.4  which  failed  to  be  classified  properly  by  the  four  com¬ 
bined  SP-LP  discriminants. 
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Plane-Wave  Array  Beam 


Adjusted-Delay  Array  Beam 


FIGURE  11-18 

SHORT  PERIOD  TIME  DOMAIN  PLOT  FOR  EVENT 

EKZ/227/01N 


•  SIN/241/  15N;  rn  =  5.  0 

b 

Surface  waves  were  not  detected  for  this  event.  The 

event  was  misclassified  by  two  of  the  six  SP  discriminants. 

The  PDE  estimate  of  5.  0  for  this  event  was  based  on  just 

one  station,  and  may  be  too  high.  The  NORSAR  m  value  for 

b 

the  event  was  4.  0. 

•  KAM/ 025/  ION;  m  =  4.6  ;  KAM/004/10N,  m  =  4.4 

b  b 

These  two  events  also  had  no  surface  wave  detection  at 
NORSAR.  However,  they  were  properly  classified  by  the  SP 
discriminants. 


YUG/ 180/  01N,  m  =4.6;  M  (Rayleigh)  =2.6,  M  (Love)  =  3.  1 
b  s  s 

This  event  was  misclassified  by  the  Mg(Rayleigh)-mb 

discriminant.  We  feel  that  this  is  because  the  m  value  re- 

b 

ported  by  PDE  may  have  been  too  high  for  this  event.  Our 

value  of  4.  6  is  the  average  of  the  four  teleseismic  PDE  stations 

that  reported  a  magnitude.  These  four  m  values  were  4.  0, 

b 

4.  0,  4.1,  and  6.  1,  respectively.  LASA  reported  m  =  4.  0  for 

b 

this  event,  and  the  NORSAR  value  computed  by  Texas  In¬ 
struments  was  3.  5.  A  value  of  mb  around  4.  0  would  imply 
correct  classification  of  this  event  by  all  discriminants. 


It  should  in  this  context  be  noted  that  we  have  had  several  simi¬ 
lar  classification  problems  with  events  for  whicft  the  PDE  mb  value  apparently 
has  been  too  high,  mainly  because  near  regional  stations  have  been  used  for  the 

PDE  determination  (Special  Report  No.  12,  1973).  In  fact,  unreliable  m  es- 

b 

timates  appear  to  be  a  major  obstacle  in  obtaining  a  satisfactory  M  -  m  sep- 

s  b 

aration  for  NORSAR  data. 


Table  II—  3  lists  the  computed  probability  of  correct  classifica¬ 
tion  for  each  individual  discriminant  based  on  all  events  of  m  >  4.4.  The 

b  ~ 
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PROBABILITIES  OF  CORRECT  IDENTIFICATION  (SHALLOW  EARTHQUAKE 
VERSUS  PRESUMED  EXPLOSION)  BY  INDIVIDUAL 


method  described  in  the  preceding  subsection  was  used  in  this  computation. 
For  the  purpose  of  comparing  the  individual  discriminants,  we  also  included 
in  Table  II- 3  the  probabilities  of  correct  classification  based  upon  only  those 
events  for  which  all  nine  discriminants  had  been  evaluated.  The  following 
main  points  summarize  the  performance  of  the  various  discriminants: 

•  The  four  combined  SP-LP  criteria  produce  a  much  better  sep¬ 
aration  between  earthquakes  and  presumed  explosions  than  the 
short  period  criteria. 

•  Probabilities  of  correct  classification  are  98-99  percent  for 
the  four  SP-LP  discriminants,  while  being  around  90  percent 
for  SP  discriminants. 

•  The  Love  wave  discriminants  appear  to  give  slightly  better  sep 
aration  than  the  Rayleigh  wave  discriminants  for  the  common 
event  set,  although  the  difference  is  only  marginal. 

•  Of  the  SP  discriminants,  the  ones  based  on  complexity  give  the 
best  separation.  However ,  the  presumed  explosion  population 
is  heavily  weighted  toward  Eastern  Kazakh  events  (distance  = 
38  degrees),  and  results  from  Special  Report  No.  11  (1973) 
show  that  spectral  discriminants  in  general  are  more  effective 
when  more  close-in  presumed  explosions  are  included. 

An  important  consideration  when  evaluating  these  discriminants 
is  of  course  how  often  they  can  be  applied.  It  appears  that  surface  waves  for 
Eurasian  presumed  explosions  cannot  in  general  be  detected  at  NORSAR  for 
events  of  mfa  less  than  5.  0.  Also,  it  seems  to  be  possible  to  detect  Rayleigh 
waves  more  often  than  Love  waves,  especially  for  presumed  explosions.  We 
have  not  tried  to  use  the  absence  of  visible  surface  waves  as  "negative  evidenc 
of  a  presumed  explosion,  although  this  is  certainly  a  way  to  extend  the  useful¬ 
ness  of  the  combined  SP-LP  discriminants. 


In  contrast  to  the  above  situation,  short  period  discriminants 
may  be  applied  at  considerably  lower  magnitudes  than  those  criteria  that  in¬ 
volve  surface  waves.  However,  our  limited  data  indicate  that  SP  discrimin¬ 
ants  do  not  perform  very  well  for  events  of  bodywave  magnitude  significantly 
below  5.  0  (Special  Report  No.  11,  1973). 

C.  PERFORMANCE  OF  MULTIPLE  DISCRIMINANTS 

In  order  to  evaluate  what  improvements  may  be  expected  from 
applying  multiple  discriminants,  it  is  important  to  know  their  interdependence. 
Obviously,  the  more  independent  two  discriminants  are,  the  more  improve¬ 
ment  can  be  achieved  by  combining  the  two. 

Table  II-4  presents  the  results  obtained  by  correlating  discri¬ 
minant  values  for  the  earthquakes  of  m  at  least  4.4  in  the  data  base.  The 

b 

discriminant  values  from  Table  II-2  were  used  for  this  computation  in  order 
to  minimize  bias  caused  by  m^  dependent  discriminants.  It  was  decided  that 
the  presumed  explosion  population  (12  events)  was  too  small  to  give  any  mean¬ 
ingful  correlation  estimates  for  this  class  of  events.  Although  it  should  be 
stressed  that  Table  II-4  only  is  intended  to  give  a  general  idea  of  the  relation¬ 
ship  between  discriminants,  the  following  points  may  be  made: 

•  Correlation  is  high  (0.  5  to  0.  9)  within  each  of  the  three  subclasses 

SP  criteria  based  on  complexity 
SP  criteria  based  on  spectral  contents 
Combined  SP-LP  criteria 

•  Correlation  is  fairly  high  (0.  3  to  0.  7)  across  the  two  classes  of 
SP  criteria 

•  Correlation  is  low  (0  -  0.  5)  between  SP  discriminants  and  com¬ 
bined  SP-LP  discriminants.  However,  there  does  appear  to  be 
a  consistently  positive  correlation  even  in  this  case. 


It  was  indicated  in  Subsection  II- A  how  the  performance  of  two 
combined  discriminants  could  be  measured.  Table  II- 5  presents  the  results 
achieved  by  combining  the  M  (Rayleigh)-  m  discriminant  with  each  of  the 

S  D 

other  criteria.  The  table  shows  that  the  probability  of  correct  identification 
typically  improves  from  about  0.  98  to  about  0.  99.  Again,  in  view  of  the 
limited  common  event  set  (only  12  presumed  explosions)  these  numbers  must 
be  interpreted  merely  as  an  indication  of  the  expected  improvements. 

It  would  be  desirable  to  present  plots  of  combined  discriminants 

in  order  to  get  a  visual  indication  of  the  improvements  achieved.  However, 

this  is  not  a  straightforward  task;  for  example,  a  plot  of  Spectral  Ratio  versus 

M  would  produce  severe  overlap  between  the  presumed  explosion  and  earth- 
s 

quake  populations.  This  is  because  both  discriminants  are  inherently  depen¬ 
dent  upon  m  ,  and  will  not  produce  good  separation  if  m  is  ignored.  Abetter 
b  b 

result  is  achieved  using  the  AL  and  AR  discriminants,  which  have  already 
been  normalized  with  respect  to  m^.  Figures  11-19  through  11-22  present 
combinations  of  these  two  discriminants  and  two  SP  criteria  (Log  P30  and 
Spectral  Ratio). 

These  plots  are  somewhat  inconclusive  with  respect  to  the  po¬ 
tential  of  combined  discriminants,  since  they  do  not  demonstrate  any  definite 

improvement  in  separation  compared  to  the  individual  plots  for  AR/m,  and 

b 

AL/ m^  (Figures  II-4  and  II - 5 ) .  It  appears  that  the  main  contribution  from 
SP  discriminants  would  be  to  give  evidence  to  corroberate  indications  of 
earthquake  or  presumed  explosion  already  provided  by  the  SP-LP  criteria. 

In  addition,  of  course,  the  SP  discriminants  are  valuable  in  those  cases  when 
surface  waves  cannot  be  detected. 


TABLE  II- 5 


PROBABILITIES  OF  CORRECT  IDENTIFICATION  (SHALLOW  EARTHQUAKE 
VERSUS  PRESUMED  EXPLOSION)  FOR  DISCRIMINANTS  COMBINED 

WITH  M  (RAYLEIGH)  -  mL 

S  D 


Discriminant 

Individual 

Discriminant 

Combined  with 

M  (Rayleigh)-  m 

S  0 

1 

Log  P30 

0.  92 

0.  985 

2 

Log  Autocorrelation 

0.  92 

0.  99 

3 

Log  Envelope  Difference 

0.  93 

0.  985 

4 

Dominant  Period 

0.  85 

0.  99 

5 

Spectral  Ratio  3/  1 

0.  92 

0.  99 

6 

Spectral  Ratio  3/2 

0.  90 

0.  99 

7 

M^Rayleigh)  -  m^ 

0.  98 

0.  98 

8 

M  (Love)  -  m, 
s  b 

0.  99 

0.  99 

9 

AR/m, 

b 

0.  985 

0.  985 

10 

AL/  m, 
b 

0.  99 

0.  99 

Log  AR  (Rayleigh) 
FIGURE  11-19 

AR  VERSUS  P30  DISCRIMINANT 
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Log  AR  (Rayleigh) 

FIGURE  11-20 

AR  VERSUS  SPECTRAL  RATIO  (1.  5  -  5)/(0  -  0.  55)  Hz 
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FIGURE  11-21 


AL  VERSUS  P30  DISCRIMINANT 
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Log  AL  (Love) 

FIGURE  11-22 

AL  VERSUS  SPECTRAL  RATIO  (1.  5  -  5)/(0  -  0.  55)  Hz 
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SECTION  III 


SUMMARY  AND  CONCLUSIONS 


The  capability  of  the  short  and  long  period  Norwegian  Seismic 
Array  (NORSAR)  to  discriminate  between  Eurasian  earthquakes  and  presumed 
explosions  has  been  evaluated.  The  total  data  base  consists  of  257  events,  24 
of  which  are  presumed  explosions.  All  of  the  earthquakes  are  either  shallow 
or  of  unknown  depth.  The  majority  of  the  presumed  explosions  (19)  occurred 
in  Eastern  Kazakh,  at  a  distance  from  NORSAR  of  38  degrees.  Bodywave 
magnitudes  for  the  events  in  the  data  base  ranged  from  5.  6  to  2.  7  for  the 
earthquakes  and  from  6.  3  to  4.  4  for  the  presumed  explosions. 

A  total  of  nine  standard  discriminants  have  been  analyzed.  Four 
of  these  are  combined  SP-LP  criteria,  which  compare  Love  or  Rayleigh  wave 
energy  recorded  at  NORSAR  to  bodywave  magnitude.  Three  are  SP  discrimi¬ 
nants  based  on  waveform  complexity  while  the  remaining  two  measure  the  SP 
spectral  energy  distribution. 

The  major  results  from  this  study  can  be  summarized  as 

follows: 

•  All  of  the  24  presumed  explosions  and  approximately  90  percent 
of  the  earthquakes  in  the  data  base  were  detected  on  NORSAR 
Short  Period  data. 

•  Twenty  of  the  24  presumed  explosions  and  approximately  70 
percent  of  the  earthquakes  in  the  data  base  were  detected  on 
NORSAR  Long  Period  data. 

•  Rayleigh  waves  were  detected  more  often  than  Love  waves  fot 
the  presumed  explosions  (20  detections  as  compared  to  14).  For 
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the  earthquake  population,  the  detectability  of  Rayleigh  waves 
was  only  slightly  better  than  that  of  JLove  waves. 

Each  of  the  four  discriminants  based  on  surface  wave  measure¬ 
ments  gave  markedly  better  separation  between  earthquakes 
and  presumed  explosions  than  any  of  the  short  period  criteria. 

All  of  the  four  combined  SP-LP  discriminants  produced  close 

to  complete  separation  for  our  event  set.  The  M^(Love)-  m^ 

discriminant  appeared  to  give  a  more  distinct  separation  than 

M  (Rayleigh) -m,  although  the  latter  one  could  be  applied  to 
s  b 

more  events.  The  AR/m,  and  AL/m,  discriminants  performed 

b  b 

about  equally  well  on  the  common  event  set. 

Of  the  short  period  discriminants,  complexity  criteria  worked 
slightly  better  than  the  spectral  ones.  However,  it  appears 
that  the  spectral  discriminants  perform  relatively  better  for 
close-in  presumed  explosions. 

The  combined  SP-LP  discriminants  were  not  able  to  give  posi¬ 
tive  identification  for  events  of  m  =  5.  0  and  lower,  since  sur- 

b 

face  waves  could  not  in  general  be  detected  at  NORSAR  for 
presumed  explosions  in  this  magnitude  range.  However,  if 
"negative  evidence"  (i.  e.  ,  the  absence  of  detectable  surface 
waves)  is  accepted  as  a  way  to  identify  a  presumed  explosion, 
the  NORSAR  SP-LP  discriminants  appear  to  be  effective  down 
to  approximately  mfa  =  4.  5.  This  is  because  surface  waves  are 
detected  for  virtually  all  Eurasian  earthquakes  above  this  mag¬ 
nitude. 

The  short  period  discriminants  worked  relatively  well  for  events 

of  m  close  to  5.  0.  However,  based  on  our  limited  presumed 
b 

explosion  data,  indications  are  that  our  SP  discriminants 
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cannot  be  expected  to  operate  well  at  significantly  below 
5.  0  at  teleseismic  distances. 

•  Correlation  is  high  between  discriminants  within  each  of  the 
three  subclasses  (combined  SP-LP,  SP  complexity,  SP  spectral 
content).  However,  the  combined  SP-LP  criteria  are  only 
slightly  correlated  with  the  SP  discriminants,  thus  indicating 
that  some  improvement  may  reimlt  from  multivariate  criteria. 

•  Detailed  analysis  of  some  "difficult"  events  indicates  that  there 
are  cases  where  SP  discriminants  work  better  than  one  of  the 
four  combined  SP-LP  criteria.  However,  in  view  of  the  high 
false  alarm  rate  (close  to  10  percent)  for  our  SP  discriminants, 
it  appears  that  decisions  (earthquake  versus  presumed  explosion) 
based  on  SP  information  alone  will  generally  have  a  substantial 
probability  of  error.  The  main  benefit  of  the  standard  SP  dis¬ 
criminants  thus  appears  to  be  to  provide  corroborative  evidence 
to  the  SP-LP  criteria. 


Since  all  the  discriminants  studied  in  this  report  are  m  depen- 

b 

dent,  it  is  essential  to  have  a  precise  estimate  of  this  parameter. 
It  appeared  that  in  some  cases  where  we  misclassified  an  earth¬ 
quake  the  main  problem  was  an  unreliable  m  estimate. 

b 


It  is  important  to  remember  that  these  conclusions  are  based 
on  the  performance  of  a  single  array  at  teleseismic  distances;  they  are  not 
necessarily  applicable  to  network  discrimination  performance  or  to  first-zone 
(  A  <  15°)  discrimination  capability. 
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APPENDIX  A 

DISCRIMINATION  STATISTICS 


Tables  A-l  and  A-2  list  the  discriminant  values  computed  for 
each  event  in  the  data  base  for  this  report.  These  values  are  identical  to  those 
used  for  the  plots  of  individual  discriminants  in  Figures  II-2  through  II- 11. 
Events  are  ordered  from  high  to  low  bodywave  magnitude.  Discriminant  values 
are  listed  in  the  following  order: 

1.  M  (Rayleigh  wave) 

s 

2.  M  (Love  wave) 

s 

3.  Log  AR  (Rayleigh  wave  energy) 

4.  Log  AL  (Love  wave  energy) 

5.  Log  P30  Mean  Square 

6.  Log  Autocorrelation  Mean  Square 

7.  Log  Envelope  Difference 

8.  Dominant  Period  (Seconds) 

9.  Spectral  Ratio  (1.5  -  5)/(0  -  0.  55)Hz  (Band  3/Bandl)  in  dB 

10.  Spectral  Ratio  (1.5  -  5)/(0.  55  -  1.  5)Hz  (Band  3/Band  2)  in  dB 

For  those  events  that  had  no  detection  of  Rayleigh  waves,  Love 
waves  or  short  period  P-waves,  the  corresponding  discriminant  values  have 
not  been  computed.  Also,  in  some  cases  the  occurrence  of  data  spikes  made 
it  impossible  to  apply  one  or  more  long  period  discriminants.  Finally,  it 
should  be  noted  that  M  (Love)  was  not  computed  for  any  of  the  18  earthquakes 
from  1971. 
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